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@ Prompt gamma neutron activation analysis system. 

(57) A system for determining depth profiles of 
concentrations of hazardous elements in soils 
comprises a neutron source (22) for generating 
neutrons of a first energy level and irradiating a 
volume of soil with the neutrons, where nudear 
reactions are effected within the soil and gam- 
ma radiation is emitted from the soil to an array 
of gamma detectors (24) for detecting gamma 
radiation emitted from the soil ; source elec- 
tronics (10) for controlling the width of regularly 
repeated pulses of neutrons generated by the 
neutron source ; detector electronics (1 2) as- 
sociated with the gamma detectors (24) for 
amplifying and digitaiizing signals generated by 
the gamma detectors and storing data repre- 
senting the digital ized signals ; spectral 
analysts software for analyzing the data and 
* determining the concentrations of selected 
target elements in the soil ; and acquisition 
interface electronics (14). 
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The present invention generally relates to PGNAA (prompt gamma neutron activation analysis), and more 
particularly relates to methods and apparatus that apply PGNAA to soil remediation. 

With the growing awareness of the contamination of large tracts of land with chemically or radioactively 
hazardous elements, there is a corresponding international effort to initiate remediation activities to restore 
5 affected regions to an environmental status considered acceptable. To this end, soil washing and other meth- 
ods are being developed. For these methods to be technically efficient and cost effective, it is necessary to 
accurately identify where the contamination in a field is located. In addition, it is necessary to have, to the great- 
est extent possible, inline monitoring of remediation process streams, to determine when the treated soil is 
acceptable for release back to the f ield. The latter need has been met by a system described in U.S. Patent 
10 Specification No. 5,133,901 (Peterson etal.). 

Neutron-induced reactions can be divided into two broad categories, threshold reactions and exoergic re- 
actions. For threshold reactions, energy in the form of neutron kinetic energy is required to supply a certain 
minimum energy to the reaction system before the reaction can proceed. Neutrons with energies below this 
minimum threshold energy are incapable of inducing the nuclear reaction. For exoergic reactions, the threshold 
15 is zero; that is. neutrons with all energies are capable of inducing the reaction. Since neutrons lose energy via 
nuclear collisions, the minimum energy possiblefor a neutron is determined by the thermal motion of the atoms 
in the stopping medium. Neutrons with this minimum average energy are referred to as thermal neutrons and 
have a mean energy of approximately 0.0252 eV. 

Figure 1 illustrates the process of neutron activation at a nuclear level. A neutron of energy E collides with 
20 the nucleus of an atom in the sample and initiates a reaction. For a neutron of thermal energy, the reaction 
might be absorption of the neutron into the nucleus, creating the next higher mass isotope of that element. If 
the neutron is more energetic (that is, with several mega-electronvolts of kinetic energy), other nuclear reac- 
tions are possible. These other reactions include inelastic scattering from the nucleus, exciting the atom ac- 
cording to its internal structure of quantum levels, or other reactions ((n,p). (n, alpha), (n,2n). etc.) in which 
25 nuclear transmutation to another element occurs. In each of these cases, the residual nucleus is left in a highly 
excited internal state, and decays to its ground state almost instantaneously (10~ 14 seconds or less), emitting 
a gamma ray of several mega-electronvolts of energy. 

The energy of this gamma ray is uniquely characteristic of the quantum structure of the residual nucleus, 
and thus is a signature of the original target nucleus. The number of atoms of each of the elements of interest 
30 in a sample can be estimated by detecting and collecting the spectrum of gamma rays emitted by the sample 
and integrating the appropriate peaks. 

The PGNAA process is governed by the following equation: 

A = N a <|> B, 

where: 

35 A = disintegrations per second producing the desired gamma rays, 
N = the number of target nuclei for the reaction, 
a = the reaction cross section (10- 24 cm 2 ), 
4> = the flux of neutrons of the required energy (cm-^sec- 1 ), 

B = the branching ratio, that is, a fraction between 0 and 1 indicating how often this capture produces 

40 the gamma ray of interest. 

In light of the growing need for soil remediation, a primary object of the present invention is to provide meth- 
ods and apparatus for accurately characterizing soil fields. A further object of the present invention is to provide 
PGNAA methods and apparatus for providing a depth profile of the contaminants in a soil field. Another object 
of the present invention is to provide PGNAA methods and apparatus for detecting trace elements in a soil 

45 field. 

A method for analyzing an interrogation volume in accordance with the present invention is characterized 
by the steps of: irradiating an interrogation volume with a neutron burst characterized by an intensity and pulse 
width, and thereby effecting an emission of gamma radiation from the interrogation volume; acquiring groups 
of gamma radiation data during a plurality of time intervals, each group being indicative of the number or in- 

so tensity of gamma rays and energy of the gamma rays during a corresponding time interval, and a first group 
corresponding to a first interval during which the neutron burst is on; and analyzing the gamma data to detect 
the presence and determine the amount of a prescribed target element in the interrogation volume. In one pre- 
ferred embodiment of the invention, the pulse width is approximately 100-200 us and the plurality of time in- 
tervals spans approximately 5000 ^s. In addition, in one embodiment described herein, the method includes 

55 generating neutrons with energies of approximately 14 MeV, 3 MeV, and 0.025 eV. While these energies are 
conveniently available, other energies are also available and useful for various applications, that is, 750 keV 
neutrons can be produced using a proton linear accelerator and a lithium target 

The method may also advantageously includes the step of analyzing the gamma data to determine wheth- 
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er the prescribed target element is present within a first prescribed depth range in the interrogation volume. 
Further, the method may include the step of generating depth profile data representing the amount of the target 
element at a plurality of depth ranges in the interrogation volume. 

In one example of the invention, the analyzing step for determining whether the target element is present 
5 within a first prescribed depth range comprises the step of computing a depth at which the neutrons have en- 
ergy appropriate to initiate a reaction that produces gamma radiation indicative of the target element 
In a second example of the invention, the method further comprises the following steps: 

determining calibration data by carrying out the following sub-steps: 
(1) placing a slug of mass M of the target element at a plurality of depths, including: 0 cm.; X 1f where X, 
w is the deepest depth from which gamma rays of energy can escape the interrogation volume in sufficient 

numbers to be detected; and X 2 , where X 2 is the deepest depth from which gamma rays of energy E 2 can 
escape the interrogation volume in sufficient numbers to be detected; and measuring yields Y E1 (0), Ye^(X^) 
of gamma rays of energy Ei at depths of 0 inches and X 1f respectively, and Y E2 (0), Y E 2(X 1 ) I Y^X^ of gam- 
ma rays of energy E 2 at depths of 0 cm., X, and X 2 , respectively; 
15 (2) defining the following ratios: 

performing a field measurement of yields y(Ei), y{E£ of gamma rays of energies E 1t E 2 , respectively; 

p m\ - >wa 

* 21(0) - VMO) 

n /v) - Yf^XQ 

20 " y £f <X 1 ) 

performing a depth inference calculation by determining whether y(E 2 )/y(E 1 ) is greater than R 2 i(X n ) and 
whether y(EO is greater than 0, and then determining a range of depths at which the target element is located 
by the following formulas: 

Case 1 

25 if y(E 2 )/y(E 1 ) is greater than R 21 (X,) and y(E t ) is greater than 0, then the target element is present be- 

tween 0 and Xi and between X 1 and X 2 ; 
Case 2 

if y(E 2 )/y(E 1 ) is less than R 21 (X0 and y^) is greater than 0, then the target element is between 0 and 

X,; 

30 Case 3 

if y(E£ is greater than 0 and y<E 1 ) is 0, then the target element is between X, and X 2 ; and 
then, performing an equivalent mass calculation by the following formulas: 
Case 1 



40 



ESmax(X<\ ~ *2) Y^X^ 

" mln (X 1 - X?) Yg^XO 



Case 2 



ES^fO-X,) = -£^M 

45 » Y E**V 

ES mln (D-X,1 = 



Y«(0) 



or, 



so ES max <0-X,) - -f~^M 



Y«(0) 



Case 3 



55 E S max(Xi - X2) - J^fw. M 
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(ES^ and ES^,, represent the maximum and minimum equivalent mass within the specified ranges.) 

The present invention also encompasses systems (hardware) for carrying out the steps described above. 
For example, the present invention encompasses a system for analyzing an interrogation volume, character- 
ized by: means for irradiating an interrogation volume with a neutron burst characterized yb an intensity and 

5 pulse width so as to effect an emission of gamma radiation from said interrogation volume; means for acquiring 
groups of gamma radiation data during a plurality of time intervals, each group being indicative of the number 
or intensity of gamma rays emitted and energy levels of said gamma rays during a corresponding time inter- 
val, and a first group corresponding to a first interval during which said neutron burst is on; and means for ana- 
lyzing said gamma data to detect the presence and determine the amount of a prescribed target element ins 

10 aid interrogation volume. 

The invention also encompasses a preferred mobile system embodiment, for determining depth profiles 
of concentrations of hazardous elements in soils, characterized by: a neutron source for generating neutrons 
of a first energy level and irradiating a volume of soil with the neutrons, whereby nuclear reactions are effected 
within the soil and gamma radiation is emitted from the soil; an array of gamma detectors for detecting gamma 

15 radiation emitted from the soil; source electronics means associated with the neutron source for controlling a 
pulse width of neutrons generated by the neutron source, and for controlling the interval of time between suc- 
cessive neutron pulses; detector electronics means associated with the gamma detectors for amplifying and 
digital izing signals generated by the gamma detectors and storing data representing the digrtalized signals; 
spectral analysis means for analyzing the data and determining the concentrations of selected target elements 

20 in said soil; and an acquisition interface module (AIM) for controlling the timing of the source and detector elec- 
tronics such that the neutron source generates neutrons in a burst of a prescribed pulse width and the detectors 
and detector electronics means detect gamma rays during a plurality of time intervals associated with the burst 
and acquire groups of data, each group being indicative of a number of gamma rays emitted and their energy 
levels during a corresponding time interval, and a first group corresponding to a first interval during which neu- 

25 trons are being generated by the source. A most preferred embodiment further comprises a mobile platform. 
In this embodiment, the array of gamma ray detectors comprises a plurality of gamma ray detectors mounted 
on the mobile platform. This embodiment also comprises means for determining a depth profile of a target 
element. 

In addition, the present invention encompasses a high sensitivity PGNAA method for detecting particular 
30 signature gamma rays generated by neutron- induced reactions within a sample interrogation volume. The re- 
actions include threshold or fast neutron-induced reactions, which occur instantaneously; thermal neutron-in- 
duced reactions, which require time for the neutrons to be moderated to thermal energies; and radionuclide 
production. The method is characterized by the following steps: (1 ) irradiating the sample with a pulsed neutron 
source; and (2) if the gamma rays of interest are generated by a threshold or fast neutron-induced reaction, 
35 counting gamma rays emitted from the sample only during the neutron pulse on-times, thereby enhancing de- 
tection of fast neutron-induced gamma rays relative to thermal neutron-induced gamma rays, radionuclide 
gamma rays, and background gamma rays; or (3) if the gamma rays of interest are generated by a thermal 
neutron- induced reaction, counting gamma rays emitted from the sample immediately following the neutron 
pulse on-times and for a duration of hundreds of microseconds thereafter. Other embodiments of the invention 
40 also include a step of or means for improving gamma signal-to- noise ratio. Other features of the present in- 
vention are described below. 

In order that the invention can be more clearly understood, convenient embodiments thereof will now be 
described, by way of example.with reference to the accompanying drawings, in which: 

Figure 1 is a schematic diagram of a basic PGNAA system. 
45 Figure 2A, which best illustrates the invention, is a schematic diagram of one embodiment of a PGNAA 

system in accordance with the present invention; Figure 2B is a more detailed schematic diagram of the elec- 
tronics of the system depicted in Fig. 2A. 

Figure 3 is a timing diagram for a PGNAA system in accordance with the present Invention. 

Figures 4Aand 4B depict exemplary gamma ray spectra obtained in accordance with the present invention. 
so Figure 5 is a graph of predicted number of reactions per atom per neutron versus depth. 

Figure 6 is a graph of experimentally measured and predicted fission responses for 236 U versus depth. 

Figure 7 depicts exemplary gamma ray spectra obtained in accordance with another embodiment of the 
present invention. 

Figure 8 is a graph of the intensity of the 2223 keV peaks of groups 1-8 (G1-G8) of Fig. 7. 
55 The system described herein employs PGNAA to determine the concentration and depth profile of hazar- 

dous elements in a field of soil. One embodiment employs a neutron generator tube in three configurations as 
a source of 14 MeV neutrons, 3 MeV neutrons, or thermal (-0.025 eV) neutrons at the soil surface. As the 
neutrons penetrate the soil, they activate the contents by both thermal capture and fast neutron induced re- 
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actions. The concentration of selected elements is determined by detecting the resulting instantaneously emit- 
ted (prompt) gamma rays, which are characterized by high and distinct energies indicative of the activated ele- 
ments. A depth distribution of the concentration of hazardous elements beneath the soil surface is derived 
through a knowledge of the depth at which the neutron field had the energy (thermal or some threshold) re- 
5 quired to initiate the reaction. Two methods for deriving a depth prof ile of contaminants are described below. 

One embodiment of the present invention is a mobile system for acquiring depth profiles of concentrations 
of hazardous materials in soil. Another embodiment of the present invention is a high sensitivity PGNAA system 
for detecting trace elements. These embodiments are described below. 

10 I. MOBILE SYSTEM FOR ACQUIRING DEPTH PROFILES OF CONCENTRATIONS OF HAZARDOUS MA- 
TERIALS IN SOIL 

Figures 2A and 2B depict one embodiment of a mobile PGNAA field survey unit in accordance with the 
present invention. 

15 Referring to Figure 2A, an exemplary embodiment of a mobile PGNAA field survey unit in accordance with 

the present invention comprises a source electronics block 10; a detector electronics block 12; an acquisition 
interface electronics block 14; a Compton Suppression block 16; a computer work station 18; a mobile or mov- 
able platform 20; a pair of neutron sources 22 carried by the platform 20 and capable of producing approxi- 
mately 10 7 thermal neutrons per cm 2 per second in a layer from 3 inches to 9 inches below the surface; six 

20 gamma ray detectors 24 carried by the platform 20; a display 26 for displaying measured gamma ray spectra; 
and memory 28 for memorizing calculated concentration and depth data. In the embodiment depicted in Figure 
2A f the mobile platform 20 is approximately 6 1.89 m. long by 0.91 m. wide and is carried approximately 15 
cm. above ground level. These dimensions are only examples. The interconnections between the source elec- 
tronics and the neutron sources and the neutron detectors have been omitted from Figure 2A to make the draw- 

25 ing more legible. Further details of the system of Figure 2A are shown in Figure 2B. 

Referring now to Figure 2B, one embodiment of the invention includes the following specific components: 

Pulse Processing Electronics 12 

30 High-speed, high throughput count rate preamplifier 12-1 , for example, a Canberra Model 2101 N Transis- 

tor Reset Preamplifier. 

Detector High Voltage Supply 12-2 capable of 0-3000 volts, for example, a Canberra Model 3106D High 
Voltage Power Supply. 

High count rate, high speed spectroscopy amplifier 12-3, for example, a Canberra Model 2024 Spectro- 
35 scopy Amplifier. 

High Speed ADC 12-4 (either multiple approximation or high frequency (at least 400MHz) Wilkinson-type 
ADC), for example, a Canberra Model 8077 ADC. (A gate signal line 12-5, AIM module input 12-6, ADC input 
12-7, and PUR/LTC (Pileup Correction/Live Time Correction) input 12-8 are also depicted.) 

40 Acquisition Interface Electronics 14 

Special Timing Module 14-1 to coordinate firing of neutron generator with gamma pulse acquisition. 
Acquisition Interface Module 14-2, for example, Canberra Model 556. 

45 Compton Suppression Detector and Associated Electronics 16 

» 

23 cm diameter x 23 cm high Nal(T1) crystal with an 8 cm diameter annulus and associated photomultiplier 
tubes (see 16-1), for example, a Solon Technologies Special MBW Assembly. 

High voltage power supply 16-2, for example, a Canberra Model 3002D High Voltage Power Supply. 
so High speed amplifier 16-3, for example, a Canberra Model 2111 Amplifier. 

Constant Fraction Discriminator 16-4, for example, Canberra Model 2126. 
Gate Pulse Delay Unit 16-5, for example, Canberra Model 2055. 
Coincidence Analyzer 16-6, for example, Canberra Model 2040. 

55 Multi-Channel Analyzer and High-speed Minicomputer (block 18) 

For example, Canberra VAX/VMS Genie Workstation. 
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Neutron Generator 22 

Pulsed D-T or D-D electronic neutron generator with variable pulse widths in the 10f>750us range and 
capable of delivering 50-1000 pulses per second. For example, a specially modified MF Physics A-320 neutron 
5 generator. 

Gamma Detector 24 

20% N-type neutron resistant high purity germanium (HPGe) detector, for example, a Canberra Model 
10 GR201 9 HPGe Detector. 

7.5 liter LN 2 Dewar, for example, a Canberra 7906- 7.5SL vertical integral cryostat. 
The development of the system depicted in Figures 2A and 2B involved determining an acceptable con- 
figuration of neutron sources and gamma ray detectors, and determining the calibration of peak intensity to 
elemental concentration for each element to be detected. The system components were selected according to 
15 the elements to be detected and the reactions to be used in detecting them. For example, activation through 
thermal neutron capture is useful for detecting uranium, thorium, cadmium, copper and many other elements. 

To achieve a depth profile of an element detectable by thermal neutron activation, one should use several 
sources which produce significant thermal neutron fluxes at very different average depths in the soil. To this 
end, the system of Figures 2A t 2B scans the ground using several different types of source, including: 
20 1 ) a DT neutron generator tube, which accelerates deuterons into tritium to make d+t fusion, producing neu- 

trons of energy approximately 14 MeV; this produces a significant thermal neutron flux at approximately 
30 to 60 cm. below the soil surface; 

2) a DD neutron generator tube, which uses d+d fusion to produce neutrons of energy approximately 3 MeV; 
this produces a significant thermal neutron flux at about 5 cm. below the soil surface; 
25 3) a DD or DT neutron generator tube, surrounded above ground by a large quantity of hydrogenous ma- 

terial, which has the effect of creating a large thermal neutron flux at the soil surface; 
4) other systems (for example, the proton linear accelerator approach) capable of providing a controllable, 
regular pattern of neutron bursts. 

In the DT tube case, there will be a significant neutron flux at energies (several MeV and higher) above 

30 the respective thresholds for inducing fast neutron reactions. Although fast neutron reactions are characterized 
by cross sections considerably smaller than thermal neutron reaction cross sections, fast neutron reactions 
also produce prompt signature gamma rays of high energy and thus can add information about the depth profile 
of particular elements. This is discussed further below. 

To estimate the performance of the system depicted in Figures 2A, 2B, a case (achievable with the DT 

35 or moderated DT configuration mentioned above) where there is a region of thermal flux of approximately 10 7 
n/cm^-sec in a volume beginning 7.6 cm. below the surface and extending to 22.9 cm. below the surface was 
considered. This is achievable by using two or three neutron tubes each generating between 1 0 9 and 1 0 11 neu- 
trons per second, (it should be noted that, in many applications, lower neutron fluxes will provide adequate 
system performance, as indicated by laboratory measurements. This alleviates the level of throughput-related 

40 complexity that could be associated with a very high thermal flux.) Table 1 is based on the 1 .89 m. x 0.91 m., 
six detector embodiment depicted in Figure 2A; it lists the total number of counts (gamma rays) per second 
expected to be intercepted by the set of 6 detectors for the elements listed. In estimating the count rates listed 
in Table 1, a specific group of activation gamma rays was assumed detected for each element. The assumed 
concentrations of elements in the soil (ppm = parts per million) are consistent with guidelines for clean soil 

45 currently recognized by regulatory agencies. In addition, by assuming that, on average, background counts 
are four times the true signal, the respective measurement times required for a standard deviation of 25% in 
the inferred concentrations of the given elements has been estimated. This corresponds to 144 true counts 
recorded during the measurement interval. Further, it has been assumed that Pb will be detected through the 
use of fast neutron induced gamma rays from inelastic scattering. 

50 
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TABLE 1 



10 



15 



Element 


ppm in soil by mass 


Expected counts/sec 


Counting time (sec) 


CI 


300 


86.4 


2 


Cu 


250 


4 


36 


Cd 


10 


5.1 


29 


Fe 


300 


5.3 


28 


Hg 


2 


2 


72 


Ni 


300 


7.2 


20 


Pb 


50 


0.9 


160 


Th 


105 


3.3 


44 


U 


42 


0.8 


180 



20 



25 



30 



35 



40 



45 



Table 1 indicates that the present invention provides a substantial amount of data in minutes of interroga- 
tion time, even considering only the thermal neutron activation component of the available data. Thus, with 
an online computer analysis system, an entire 15.2 m. x 15.2 m. area can be evaluated at a single depth in 
approximately ten hours. A three region depth profile can be acquired in approximately thirty hours. These 
performance estimates are for the 1 .89 m. x 0.91 m., six detector embodiment of Figure 2A. A larger area of 
coverage per interrogation could be achieved simply by increasing the dimensions, number of detectors, and 
number of sources, in addition, it should be noted that, although the above data was generated for a few ele- 
ments of interest, the present invention may be employed to detect concentrations of elements across the per- 
iodic table. 

Both fast neutron induced reactions and thermal neutron induced reactions are mentioned above but, in 
the estimate of performance, only thermal neutron capture is considered, except in the case of Pb, for which 
inelastic scattering was assumed. As discussed above, even in nuclides with low thermal capture cross sec- 
tions (for example, 12 C or the Pb isotopes), inelastic scattering of a neutron will produce a prompt signature 
gamma ray and thus provide an alternate method of element identification. In many nuclides, threshold reac- 
tions such as (n,p) or (n,a), which often require a minimum incident neutron energy before the reaction can 
proceed, are also likely. Each of these reactions produces its own set of unique, signature prompt gamma rays. 
Therefore, extension of the neutron field to contain fast, thermal and intermediate energy components will ex- 
tend the versatility of the system by extending the variety of reaction channels that can be used for element 
identification. 

Moreover, since depth profiles of the various energy components of the neutron field can be reliably cal- 
culated with several well-validated transport codes (for example, DOT or MCNP), the response of signature 
gamma rays from reaction channels initiated by neutrons of differing energy regimes can also be used to Infer 
contamination depth profiles. This is discussed further below. 

To generate neutrons of energies sufficient to initiate most fast neutron induced reactions, the most com- 
monly used sources are electronic neutron generator tubes of the type developed for oil well logging applica- 
tions. These tubes contain a compact accelerator which propels deuterons into tritium or deuterium at an inci- 
dent kinetic energy of 100 - 200 keV. This produces fusion reactions, with a resulting neutron yield. The D + T 
reactions produce 14 MeV neutrons and the D + D reactions produce 2-3 MeV neutrons. 
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II. TIMING OF NEUTRON PULSES AND GAMMA ACQUISITION 

Whether using a Cf (Californium) source or an electronic neutron generator, there is a tendency for the 
signal processing circuitry to experience saturation effects due to the high rate of prompt gamma ray production 
while the neutrons are being emitted. This can lead to a severe reduction of gamma energy data throughput, 
to the point of rendering the acquisition system completely inactive in extreme cases. This situation can be 
remedied in two ways: (1) by employing fast electronics, including, for example, transistor reset preamplifiers 
and fast spectroscopy amplifiers having the best possible signal handling capability; and (2) by modifying the 
timing of the neutron emissions from the tube such that the rate of emissions will be matched to the circuit 
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signal processing capability. 

Most neutron generator tubes produce a very intense burst of neutrons within the first few microseconds 
of each pulse period. Atypical pulse duration is 8 \xs. This timing is not within the user's control in most systems. 
By allowing the neutron on-time (pulse width) and period to be adjustable, the same average number of neu- 

5 trons can be produced per unit time using a longer on-time per pulse period. This enables one to reduce the 
instantaneous load on the signal processing system during the neutron on-time to a level the detector elec- 
tronics can handle without severe degradation of data throughput With the present invention, the prompt gam- 
ma rays produced by fast neutrons are detectable at levels sufficient for element identification instead of being 
lost because of saturation effects in the electronics. 

10 There are two reasons why it is advantageous to separate acquisition of prompt gamma rays produced by 

fast neutrons and prompt gamma rays produced by thermalized neutrons: First, having a separation of fast 
neutron induced gamma rays and thermal neutron induced gamma rays simplifies spectral interpretation. Sec- 
ond, this optimizes the signal-to-background ratio for the thermal neutron induced prompt gamma rays by min- 
imizing the data representing long-lived activity. 

15 Such timing control requires three major system components: 

1 . A neutron generator with adjustable pulse period and neutron on-time, or duty cycle. 

2. Fast electronics in preamplifiers, spectroscopy amplifiers and other signal processing components. 

3. An electronic control system that can control the neutron pulse characteristics and also acquire and store 
the data in separate "fast neutron prompt gamma," "thermal neutron prompt gamma" and "long lived ac- 

20 tivity" buffers. 

The system must acquire gamma ray spectral data in distinct time slices: one or more representing the 
on-time of the neutrons (intentionally kept small compared to the expected time -typically several hundred mi- 
croseconds- for neutrons to thermalize (slow down) in the irradiated medium) and one or more encompassing 
some or all of the remaining time between neutron bursts. Such a scheme allows one to empirically identify 

25 the separate time regimes appropriate to the fast neutron prompt gamma rays, thermal neutron prompt gamma 
rays, and long lived activity, without including unnecessary gamma data in the respective buffers. 

This approach is illustrated in Figure 3, which illustrates the timing of the gamma ray acquisition protocol. 
As shown, the neutron burst pulse width is approximately 100 to 200 microseconds and the pulse repetition 
period is approximately 5000 microseconds. The neutron burst on-times are represented by the shaded regions 

30 in Figure 3. During the burst on-time, gamma rays produced by fast neutrons are detected; thereafter, in an 
interval encompassing approximately 1000 to 1500 microseconds starting at the end of the neutron burst on- 
time, gamma rays produced by thermal neutron activation are detected; then, after this second acquisition re- 
gime, gamma rays produced by long-lived neutrons are detected. One embodiment of a pulsed PGNAA system 
in accordance with the present invention includes the above-described neutron generator tube and detection 

35 and control electronics. The pulsed PGNAA system can advantageously be employed in the mobile field survey 
application described above. 

Data taken using this system is shown in Figures 4A and 4B. Figure 4A depicts a series of eight gamma 
ray spectra, referred to as groups G1 through G8. Each group represents a spectrum of gamma rays detected 
during an interval of 500 microseconds. Group G1 corresponds to the neutron generator on-time, that is, the 

40 first 500 microseconds. The following groups (Groups G2-G8) represent successive 500 microsecond intervals 
following the neutron generator on-period. Figure 4B depicts an expanded view of Group G2. Hydrogen and 
aluminum activation lines are identified in this figure. As mentioned above, the specific energy levels corre- 
sponding to particular elements of interest are known in advance from direct measurement or compilations of 
nuclear data. 

45 

III. INFERENCE OF CONTAMINATION DEPTH PROFILES 

Two methods for determining the depth distribution of contaminants will now be described. The first meth- 
od makes use of the known behavior of the neutron energy distribution with depth as the neutrons propagate 
so beneath the surface; the second method makes use of the known attenuation properties of soO for gamma 
rays of very different energies, emitted with known yields, from the same element 

A. METHOD BASED ON SPATIAL AND ENERGY DISTRIBUTION OF NEUTRON FIELD 

55 As mentioned above, a variety of reactions will be initiated for any given element if one employs a neutron 

field with fast, thermal, and intermediate neutrons. Since the prompt gamma rays are signatures of the target 
element and the reaction, the measured response of one of these gamma rays represents a folding of the con- 
taminant concentration profile with the spatial distribution of neutrons whose energy is appropriate to initiate 
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the reaction. Thus the depth distribution of a contaminant of interest can be determined if the spatial and energy 
distribution of the neutron field is known. In other words, since the energy of the signature gamma ray indicates 
both the element and the reaction that produced the gamma ray, the depth of the element can be bracketed 
by calculating the region of depth at which the neutron field had the energy appropriate to initiate the reaction. 
5 Since the chemical composition, bulk density, and approximate water content of soil can be measured, it 

is possible to use any one of several well known, internationally recognized transport codes (for example, DOT 
or MCNP) to calculate the neutron energy and spatial distribution of neutrons. One can also perform calcula- 
tions for a given set of conditions, and develop empirical corrections for the key parameters, for example, den- 
sity and water content. For most soils, the relative amounts of silicon, aluminum, calcium, magnesium, potas- 
10 sium, and other oxides composing most of the soil mass will not significantly affect the neutron behavior 
(though they could produce very different responses of emitted gamma rays). 

An illustration of the predictive capability of transport codes is provided by Figure 5, which shows the depth 
distribution of reactions per atom per neutron predicted with a DOT calculation for a model of a large laboratory 
mockup of soil. The code was given response functions for several reactions, including neutron induced fission 
15 of and ^ZTh, which, when folded with the energy spectrum of neutrons at each node in the model, pro- 
duced the depth distributions shown. These particular responses were selected for two reasons: 

1. They accentuate the effect of very different regimes of neutron energy; the fission cross section 
is almost entirely thermal (neutron energy below 0.5 eV) and the ^mi cross section exhibits a neutron 
energy threshold of approximately 1.3 MeV. 
20 2. The use of Solid State Track Recorders (SSTRs) containing micro deposits of 235 U enabled a direct meas- 

urement of these responses as a function of depth. 
, Figure 6 is a graph of experimentally measured ^HJ fission response and the prediction from DOT. This 
graph shows that the DOT model is capable of producing accurate predictions of the spatial distribution of neu- 
trons. In practice, it would be desirable to ensure that an appropriate placement of neutron sources is made 
25 to ensure optimized lateral uniformity of the contours of neutron flux. 

B. METHOD BASED ON DIFFERENTIAL ATTENUATION OF GAMMA RAYS 

When a nuclear reaction occurs in a sample, prompt gamma rays of several different energies are often 
30 emitted from the sample; each has a known probability of occurrence, or branching ratio (often defined as the 
number of gammas of a particular energy emitted per 1 00 reactions). This is discussed previously. Any medium 
through which the gamma rays travel will attenuate the gamma ray intensity to a degree that is a function of 
their energy. In particular, the attenuated gamma ray intensity l(X) (photo ns/cm^sec) at a depth X in an ab- 
sorbing medium, assuming a plane source of radiation of initial intensity l 0 , is given by: 
35 l(X) = l 0 e ■ 

where the attenuation coefficient n is a decreasing function of photon energy. Thus, if a nucleus emits a photon 
of energy E, with branching ratio and a photon of energy E 2 with branching ratio B 2 , the ratio of the respective 
measured yields, Y x and Y 2t after attenuation through a thickness X of soil, is given by: 

where 

Hi = \i(EJ, 
u 2 = nfEz). 

Therefore, one can determine X from the equation: 
45 ' x = Incyj/Yz) - InCBy^Bz) 

H2 - m 

Thus, the depth at which an emitting plane source is located can be inferred. By extending this treatment 
to a continuous distribution (for example, a superposition of point sources), the appropriate software can infer 
bounding concentrations and spatial distributions of a contaminant element from the intensity and ratio of meas- 
50 ured yields of gamma rays of differing energies. 

Example: Depth Inference From Yields of Two Reactions 

Consider a case where thermal reactions induce two gamma peaks of energies Ej and E 2 , respectively, 
55 where E 2 » E, (e.g., E 2 = 4-9 MeV, E, = 0-1 MeV). 
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Stop 1: Calibration 



A slug of mass M is buried at several depths, including: 

- Ocm.; 

5 - X 1f where X^ is the deepest depth from which gamma rays of energy E, can escape the soil in sufficient 

numbers to be detected; 

- X 2 , where X 2 is the deepest depth from which gamma rays of energy E 2 can escape the soil in sufficient 
numbers to be detected. Of course, X 2 is greater than X v 

Calibration measurements determine the yield Y of gammas of E n at depth X,. This yield is defined as 
w Y En (X|). Thus, Y E1 (0), Y^Xi), Y^O). Y^), Y^XJ will be known (Y^fXJ will be zero). 
The following ratios are defined: 

Ve2(0) 



R 2 ,(0) = 



d /y ) — Y E2(Xi) 

15 Y £ f(>vi) 

(Note: For E A approximately equal to 1 MeV and E 2 approximately equal to 4 MeV, X-, will be about 15 cm. and 
X 2 will be about 56 cm.) 



Step 2: Field Measurement of Gamma Ray Spectrum 

20 

The field measurement of the spectrum provides two measured yields; these are denoted y(Ez) and y(E 1 ). 
All yields wBI preferably be normalized to the number of neutrons emitted by the neutron generator that is, 
measured yields will be expressed in terms of counts per neutron. 

25 Step 3: Depth Inference Calculation 



It is important to note that the ratio y(E 2 )/y(E 1 ) will increase with depth below the surface. The location of 
the element of interest, in terms of a range of depths, is inferred by first answering the following two questions: 
lsy(E 2 )/y(Ei)<R 2 i(X 1 )? 
lsy(E 1 )<0? 

Case 1 : Element between 0 and and X^ and X 2 

y(E 2 )/y(E 1 ) wBI be greater than R 2 i(Xi) and y(E t ) will be greater than 0. 
Case 2: Element between 0 and X t 

y(E 2 )/y(E 1 ) will be less than R 2 i(Xi) and y(E 1 ) will be greater than 0. If y(E,) is less than 0, the wrong 
element has been assumed. 

Case 3: Element between X, and X 2 

y(E 2 ) will be greater than 0 and y(E0 will be 0. 



40 



Step 4: Equivalent Mass Calculation 

Case 1 : Element between 0 and X< and X A and X 2 



ES max (0-X0 = -^§-M 

45 = -g^M 

YeK*2) 

50 r E*X\) 

(Note: The product y(E 1 )R 21 (X 1 ) represents the maximum yield of E 2 gammas from 0 to X,; y(E0R 21 (0) repre- 
sents the minimum yield of E 2 gammas from 0 to X v ) 

Case 2: Element between 0 and X, 
Alternatively, the following equations can be used to calculate the maximum and minimum equivalent element 
55 slugs between 0 and 
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ES max (0-X,) = ^-W 



Case 3: Element between X-, and X 2 



ES max (Xi - = y^x^ 

ESm\n(X<\ - X2) ~ J^ v v \^ 

(ES max and ES mln represent the maximum and minimum equivalent mass within the specified ranges.) 

It should be noted that the maximum and minimum mass content can always be determined for two zones 
(0 ^ depth ^ X 1? X n ^ depth ^ X2) by using two predetermined energy levels. The analysis can be extended 
to more zones by using more energy levels. 

Next, to convert the mass content to parts per million (ppm), the calculated mass must be divided by the 
mass of the contributing soil. For example, this calculation can be performed by assuming that the shape of 
the contributing mass is a trapezoidal cylinder. The details of this calculation will be apparent or readily available 
to those of ordinary skill in the art 

Table 1, previously set forth, lists expected counts per second for a number of elements of interest. The 
table assumes typical recommended "clean" levels of the elements. In addition, a typical soil slurry (1 .76 gm/cc 
dry, 94% mass fraction powder in slurry) is assumed. Note that, for Pb, inelastic scattering of energetic neu- 
trons is assumed, with a conservative cross section assumed for the expected energy spectrum. Thermal cap- 
ture is assumed for all other elements. (This does not preclude the use of the disclosed methods and apparatus 
to measure the concentration profile of Pb using thermal capture.) Further, the time required to detect the as- 
sumed concentration with a standard deviation of 25%, amidst a background which is four times the true count 
rate, is also listed. The results listed in Table 1 demonstrate that a PGNAA system employing the neutron pulse 
and acquisition timing control described above is a powerful and practical tool. 

IV. HIGH SENSITIVITY PGNAA 



The present invention may be employed for the detection of trace elements in a sample. This capability is 
important for a number of applications, including in- srtu PGNAA analysis of soils for trace contaminants and 
on-line monitoring of soil washing (see U.S. Patent Specif ication No. 5,133,901). Several modifications of the 
methods described above have been devised to provide a trace element analysis capability. These modifica- 
tions are described below. 

Except for 4 He, all stable nuclides can capture a thermal neutron with a release of energy in the form of 
gamma rays. Since reaction cross sections generally increase with decreasing energy, thermal neutron cap- 
ture gamma rays are usually produced in the highest yields in neutron irradiations. However, background gam- 
ma rays and gamma rays from threshold reactions and activation product radionuclides will also be present 
Methods for enhancing the detection of the gamma rays of interest are required to successfully use PGNAA 
for trace element analysis. This is because the most abundant gamma rays will generally be produced by the 
most abundant elements. 

The detection of particular signature gamma rays can be facilitated by exploiting the time sequence of neu- 
tron-induced reactions. Threshold reactions (fast neutron- induced reactions) occur instantaneously; that is, the 
gamma rays are emitted within 10~ 14 seconds of neutron capture (no time is required for the neutrons to be 
moderated to lower energies). Therefore, even with the fastest available electronics, neutrons will still be pres- 
ent when the gamma rays are detected. Thermal neutron- induced reactions also proceed rapidly but time is 
required for the neutrons to be moderated to thermal energies. Typically, the time required is on the order of 
0.1-1 jas and, once formed, the thermal neutrons can have mean lifetimes of tens to hundreds of microseconds, 
depending on the neutron capture characteristics of the moderating media. During this time, thermal neutron- 
induced reactions can take place until all the neutrons are captured. Gamma rays from radionuclide production 
and background gamma rays will essentially be constant on a microsecond to millisecond time scale. Thus, if 
the nuclide of interest is to be detected via a threshold reaction, use of a pulsed neutron source and counting 
only during the neutron pulses will enhance detection of these gamma rays relative to thermal neutron-induced 
reaction gamma rays, radionuclide gamma rays, and background gamma rays. Similarly, detection of thermal 
neutron- induced reaction gamma rays can be enhanced by counting immediately following neutron bursts and 
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for a duration of hundreds of microseconds. 

High-speed detection electronics are required to exploit the time sequence of neutron-induced reactions 
and obtain increased gamma ray detection sensitivities. In a standard High-Purity Germanium (HPGe) gamma 
ray detection system, count rate throughput is controlled by: 
5 (1) charge collection times for pulses produced by gamma rays in the HPGe detector; 

(2) pulse handling characteristics of the preamplifier; 

(3) pulse handling characteristics of the amplifier; 

(4) speed of the analog-to-digital converter (ADC); and 

(5) storage time of the multichannel pulse height analyzer. 

10 The system amplifier presents the most serious bottleneck to count rate throughput However, if the HPGe 

detector is equipped with a high-speed transistor reset preamplifier and an advanced ADC (either fixed con- 
version time or Wilkinson-type) is used, count rate throughputs of up to 100,000 counts per second (cps) are 
attainable. 

At a count rate throughput of 100,000 cps, the system dead time (that is, the time required for the elec- 
ts tronics to process pulses) will be 60-80% of the clock time, depending on spectrum characteristics. Therefore, 
a pulsed neutron generator with extended neutron pulse duration is employed to decrease the neutron intensity 
per unit time and thereby enhance detection of gamma rays from fast neutron-induced reactions. The neutron 
generator produces neutron bursts with durations in the 100-700 \is range. The neutron output of this generator 
is about 5 x 10 7 neutrons/cm 2 /second. Control of the generator is accomplished through an acquisition interface 
20 module (AIM). 

The AIM is designed to run the neutron generator at a preselected pulse length and frequency while at 
the same time controlling the multichannel analyzer (MCA) to obtain gamma ray spectra for predetermined 
acquisition periods synchronized with the operation of the neutron generator. Thus, for example, the neutron 
generator may be run for 200jis while the MCA obtains a synchronized 200|is spectrum followed by a series 

25 of seven 400 us spectra. Fast neutron- induced reaction gamma rays would be detected only in the first counting 
interval whereas thermal neutron- induced reaction gamma rays would be detected mainly in the first four 
groups. Background gamma rays and gamma rays emitted by radionuclides produced by neutron activation 
would be detected with equal probability in all eight intervals. 

Figure 7 depicts an example of high-sensitivity PGNAAdata acquired with high-speed timing techniques. 

30 The neutron generator was run for 500us at a rate of 200 pulses per second. Each spectrum corresponds to 
a counting period of 250ns. The first two spectral groups are thus encompassed by the generator-on period. 

Gamma rays from fast neutron- induced reactions appear only in the first two groups (that is, the first two 
spectra), which contain many more counts than the following six groups. Both groups taken with the neutron 
generator on, correspond to very high dead times (high count rate losses). Although thermal neutron-induced 

35 reaction gamma rays are present in the first two groups, their intensity is low because of losses and the fact 
that the thermal neutron reaction rate had not yet peaked. A prominent peak in the spectra, at 2223 keV, is a 
consequence of the capture of thermal neutrons by hydrogen. Figure 8 shows a plot of the intensity of the 
2223 keV gamma rays for each of the spectral groups. The maximum count rate occurs in group 3 (the first 
250us after the neutron pulse ended), and the intensity falls off exponentially with time in the next few groups. 

40 The remaining signal in groups 7 and 8 is due to so-called "wall return" thermal neutrons. (Wall return is a lab- 
oratory effect.) The intensity of the 2223 keV gamma rays decreases by a factor of 2 every 150us, a rate that 
is determined by the neutron capture characteristics of the thermalizing medium, which in this case was soil 
below. the neutron generator. Groups 7 and 8 contain mainly "long-lived" gamma rays from background and 
from nuclides produced in neutron activation (such as 16 N with a half life of 7 seconds produced by the 

45 16 0(n,p) 16 N reaction on oxygen in the soil). These gamma rays are present at equal intensity in all of the groups. 

It can be seen that use of the data from groups 1 and 2 will result in much better sensitivity for fast neutron- 
induced reactions because only the background that occurs in the narrow time interval corresponding to the 
generator-on period is encountered. Similarly, maximum sensitivity for neutron capture gamma rays results 
from summing the data in groups 3 and 4. 

50 The methods described herein offer a number of improvements to basic PGNAA techniques, including: 

(1) increased sensitivity due to higher count rate input; 

(2) increased sensitivity due to detection of reaction gamma rays during the time interval in which they 
occur (that is, increased signal-to- background); and 

(3) separation of gamma ray events from different types of reactions. Fast neutron- induced reactions and 
55 thermal neutron- induced reactions will occur with different spatial distributions in the sample. Forexample, 

if 14 MeV neutrons are being attenuated by a solid such as soil, the fast neutron flux will drop off expo- 
nentially with distance from the surface of the sample, whereas the thermal neutron flux will increase and 
reach a maximum roughly six to ten inches from the surface. Thus, different portions of a non homogeneous 
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sample can be sampled simultaneously. 

The sensitivity of the method may be further enhanced by the addition of an anti-Compton suppression 
ring of detectors around the HPGe detector(s). These improvements to PGNAA make possible the use of the 
technique for trace element analysis. 

5 

V. PGNAA SYSTEM FOR ANALYSIS OF MAJOR AND TRACE COMPONENTS OF OBJECTS, PACK- 
AGES OR FLOW STREAMS 

The methods discussed above are useful for analyzing liquids, soil slurries, powders, solid objects, soil 
10 cores, and packages of various sizes. In addition, the system is transportable, so that, if needed, it could be 
assembled onsite at a remote facility to characterize suspect samples without shipping the samples offsite. 
Further, as discussed above, the system is capable of real-time measurement and providing results within min- 
utes. 

Assume, for example, that the system has an interrogation volume of one cubic foot and that, using a 14 
15 MeV tube source, the moderating material produces within the interrogation volume a thermal neutron flux of 
10 7 neutrons/cm^sec and a flux of E > 1 MeV neutrons of 10 7 neutrons/cmZ-sec. Further, assume that there 
are two high purity germanium detectors respectively placed on opposite faces of the interrogation volume. 
The following table lists expected counts per second from each listed element, at typical recommended "clean" 
levels in soil. A typical sofl slurry (1.76 gm/cc dry, 94% mass fraction powder in slurry) is assumed for this 
20 illustration. Note that, for Pb, inelastic scattering of energetic neutrons is assumed, with a conservative cross 
section assumed for the expected energy spectrum. Thermal capture is assumed for all other elements (more 
recently obtained data indicate that thermal neutron capture for Pb may lead to tower detection limits for this 
element). In addition, the time required to detect the assumed concentration with a standard deviation of 25%, 
amidst a background which is four times the true count rate, is also listed. 

25 

TABLE 2 
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Element 


ppm in soil by mass 


Expected counts/sec 


Counting time (sec) 


CI 


300 


18.3 


7.9 


Cu 


250 


0.8 


180 


Cd 


10 


1.1 


131 


Fe 


300 


1.1 


131 


Hg 


2 


0.4 


360 


Ni 


300 


1.5 


96 


Pb 


50 


0.2 


720 


Th 


105 


0.7 


206 


U 


42 


0.2 


720 



45 * Table 2 demonstrates that a system based on PGNAA with the neutron pulse and acquisition timing control 
described above can be used as a practical tool for the measurement of concentrations of hazardous elements 
in a variety of samples and geometries. Such a system will also use neutron doses which are so low (typically 
less than 10 12 neutrons) that post-irradiation residual radioactivity will not be a problem. 

In summary, PGNAA can be used as a real time method for determining the depth distribution of hazardous 

so elements in a field of soil. Using successive measurements over a field with three different neutron generator 
tube configurations (DT, DD, and DD with above ground t normalization), one can achieve a depth distribution 
due to the three respective average depths of neutron thermalization these sources provide. In addition, using 
a DT source of 14 MeV neutrons, there will be sufficient neutron energy to excite inelastic scattering and other 
"fast neutron" reactions, which also produce distinct signature gamma rays of high energy. These additional 

55 gamma rays can enhance the depth profile information, according to the region of soil depth at which the neu- 
tron average energy is above the neutron energy thresholds of the respective reactions. Utilizing measured 
ratios of gammas of different energies produced by the same nuclide further enhances the depth profile in- 
formation because of the different attenuation of intensity experienced by the different energy gamma rays. 
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A three region depth profile of a 15.24 m. x 15.24 m. area can be produced in approximately thirty hours of 
interrogation time using a system including two neutron tube sources, six germanium detectors covering a 1 .89 
m. x 0.91 m. ground surface, and an online computer analysis system. Larger systems will provide commen- 
surately lower survey times. The systems described above may also be employed to detect trace elements. 

5 

Claims 

1. A method for analyzing an interrogation volume, characterized by the steps of: 
10 (a) irradiating an interrogation volume with a neutron burst characterized by an intensity and pulse 

width, and thereby effecting an emission of gamma radiation from said interrogation volume; 

(b) acquiring groups of gamma radiation data during a plurality of time intervals, each group being in- 
dicative of the number or intensity of gamma rays and energy levels of said gamma rays during a cor- 
responding time interval, and a first group corresponding to a first interval during which said neutron 

is burst is on; and 

(c) analyzing said gamma data to detect the presence and determine the amount of a prescribed target 
element in said interrogation volume. 



20 



25 



2. The method recited in claim 1, characterized in that said pulse width is approximately 100-200 us. 

3. The method recited in claim 1 , characterized in that said plurality of time intervals spans approximately 
5000 us. 

4. The method recited in claim 1 , characterized in that step (a) comprises generating neutrons with energies 
of approximately 14 MeV. 

5. The method recited in claim 1 , characterized in that step (a) comprises generating neutrons with energies 
of approximately 3 MeV. 

6. The method recited in claim 1 , characterized in that step (a) comprises generating neutrons with energies 
30 of approximately 0.025 eV. 

7. The method recited in claim 1 , characterized in that step (a) comprises generating neutrons with energies 
of approximately 750 keV. 

35 8. The method recited in daim 1, further characterized by the step of analyzing said gamma data to deter- 
mine whether said prescribed target element is present within a first prescribed depth range in said in- 
terrogation volume. 
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9. The method recited in claim 8, further characterized in that said analyzing step for determining whether 
said target element is present within a first prescribed depth range comprises the step of computing a 
depth at which said neutrons have energy appropriate to initiate a reaction that produces gamma radiation 
indicative of said target element. 

10. The method recited in claim 8, further characterized by the step of generating depth profile data repre- 
senting the amount of said target element at a plurality of depth ranges in said interrogation volume. 



11. A system for analyzing an interrogation volume, characterized by: 

(a) means for irradiating an interrogation volume with a neutron burst characterized by an intensity and 
pulse width so as to effect an emission of gamma radiation from said interrogation volume; 

(b) means for acquiring groups of gamma radiation data during a plurality of time intervals, each group 
50 being indicative of the number or intensity of gamma rays emitted and energy levels of said gamma 

rays during a corresponding time interval, and a first group corresponding to a first interval during which 
said neutron burst is on; and 

(c) means for analyzing said gamma data to detect the presence and determine the amount of a pre- 
scribed target element in said interrogation volume. 



55 



12. The system recited in claim 11, further characterized by means for generating a neutron burst by a pulse 
width of approximately 100-200 us, and where said plurality of time intervals spans approximately 5000 
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13. The system recited in claim 11 , further characterized by means for analyzing said gamma data to deter- 
mine whether said prescribed target element is present within a first prescribed depth range in said in- 
terrogation volume. 

14. The system recited in claim 13, characterized in that said means for analyzing said gamma data to de- 
termine whether said target element is present within a first prescribed depth range comprises means 
for computing a depth at which said neutrons have energy appropriate to initiate a reaction that produces 
gamma radiation indicative of said target element. 

15. The system recited in claim 13. further characterized by means for generating depth profile data repre- 
senting the amount of said target element at a plurality of depth ranges in said interrogation volume. 

16. A mobile system for determining depth profiles of concentrations of hazardous elements in soils, charac- 



(a) a neutron source for generating neutrons of a first energy level and irradiating a volume of soil with 
said neutrons, whereby nuclear reactions are effected within said soil and gamma radiation is emitted 
from said soil; 

(b) an array of gamma detectors for detecting gamma radiation emitted from said soil; 

20 (c) source electronics means associated with said neutron source for controlling a pulse width of neu- 

trons generated by said neutron source, and for controlling the interval of time between successive neu- 
tron pulses; 

(d) detector electronics means associated with said gamma detectors for amplifying and digitalizing sig- 
nals generated by said gamma detectors and storing data representing the digitalized signals; 
25 (e) spectral analysis means for analyzing said data and determining the concentrations of selected tar- 

get elements in said soil; and 

(f) an acquisition interface module (AIM) for controlling the timing of said source and detector electronics 
such that said neutron source generates neutrons in a burst of a prescribed pulse width and said de- 
tectors and detector electronics means detect gamma rays during a plurality of time intervals associated 
30 with said burst and acquire groups of data, each group being indicative of a number of gamma rays 

emitted and their energies during a corresponding time interval, and a first group corresponding to a 
first interval during which neutrons are generated by said source. 

17. The mobile system recited in claim 1 6, further characterized by a mobile, platform, wherein said array of 
35 gamma ray detectors comprises a plurality of gamma ray detectors mounted on said mobile platform. 

18. The system recited in claim 17, further characterized by means for determining a depth profile of a target 
element 

4Q 19. A high sensitivity prompt gamma neutron activation analysis (PGNAA) method for detecting particular sig- 
nature gamma rays generated by neutron-induced reactions within a sample interrogation volume, char- 
acterized by the steps of: 

irradiating the sample with a pulsed neutron source; and 

if the gamma rays of interest are generated by a threshold or fast neutron-induced reaction, count- 
45 ing gamma rays emitted from the sample only during the neutron pulse on-times, thereby enhancing de- 

tection of fast neutron- induced gamma rays relative to thermal neutron-induced gamma rays, radionuclide 
gamma rays, and background gamma rays; 

if the gamma rays of interest are generated by a thermal neutron-induced reaction, counting gam- 
ma rays emitted from the sample immediately following the neutron pulse on-times and for a duration of 
so hundreds of microseconds thereafter. 

20. The method recited in claim 19, further characterized by the step of analyzing the collected gamma data 
to determine whether a prescribed target element is present within a first prescribed depth range in the 
sample interrogation volume, and where said analyzing step for determining whether said target element 
55 is present within a first prescribed depth range comprises the step of computing a depth at which said 

neutrons have energy appropriate to initiate a reaction that produces gamma radiation indicative of said 
target element. 
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